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Abstract 

Using the Hilbert-Schmidt theorem, we reformulate the R-matrix theory 
in terms of a uniformly and absolutely convergent expansion. Term by term 
differentiation is possible with this expansion in the neighborhood of the 
surface. Methods for improving the convergence are discussed when the R- 
function series is truncated for practical applications. 

I. Introduction 

Since 1947, boundary condition methods (BCM) have played an impor- 
tant role for many quantum mechanical problems [1-37]. In the BCM formu- 
lation, configuration space is divided into two parts: internal and external 
regions. In the external region, the interaction is usually known and in many 
cases the effective two-body equation is exactly solvable. A boundary con- 
dition matrix is defined in terms of the independent external wave functions 
and their derivates at a boundary. From this information (boundary con- 
dition matrix) and the known solution in the external region, the S-matrix 
and the cross-section can be calculated. There are two boundary condition 
matrices: R matrix and P matrix. The R-matrix, which is the inverse of the 
logarithmic derivative of the external channel wave function at the surface, 
was first introduced by Wigner and Eisenbud in 1947 [1]. A detailed account 
of the R-matrix theory of nuclear reactions is given in [3] . The P-matrix is the 
inverse of the R-matrix. The P-matrix formulation of nuclear reactions has 
not been used extensively except for the nucleon-nucleon scattering problem 
[5-10]. 

The R-matrix theory is extensively employed for describing energy de- 
pendence of the cross-section for various binary nuclear processes including 
both elementary and composite nuclear particles [3, 11-17] and is often used 
to extrapolate experimental data for the cross-section which are dominated 
by the contributions from a few resonance or bound states. 
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In atomic physics, the exchange interaction, which is the most difficult 
part of the problem to calculate accurately, is only in the internal region, 
and interactions in the external region reduce to long range local poten- 
tials [18]. The R-matrix code is a very powerful computation code [19, 20] 
for calculating electron-atom collisions and photoabsorption processes. For 
example, the inner-shell photoionization cross-sections calculated by the R- 
matrix code [21, 22] are in excellent agreement with the recent experimental 
measurements [38]. 

The R-matrix method for studying low-energy electron- molecule collisions 
was developed in [23-25]. It has been used to describe elastic scattering, 
electronic excitation, vibrational excitation [27] and dissociative attachment 
[29]. 

The R-matrix method is based upon expanding the total wave function 
\& for any energy in the internal region in terms of the complete set of eigen- 
f unction Xx of that region, 

* = ^ AxXx, (1) 

A 

where Xx are defined by the equation 

HXx = ExXx, (2) 

and satisfy certain R-matrix boundary conditions on the surface [3]. It is 
known that there is a formal problem with the wave function expansion used 
in the conventional approach for the R-matrix theory. Either the expansion 
given by Eq. (1) is not uniformly convergent in the neighborhood of the 
surface, or term-by-term differentiation of the expansion is not admissible 
[3,18]. To avoid these difficulties variational formalisms [30-36] were proposed 
with basis functions which do not satisfy the R-matrix boundary conditions. 

The R-matrix theory is rigorous and therefore, there remains, from the 
formal point of view, a problem of how to formulate the R-matrix method in 
terms of a uniformly and absolutely convergent expansion. 

In this paper we present solutions of this formal problem and discuss 
various approximations of the R-function. In Section II, we describe in some 
detail the conventional formulation of R-function in terms of the expansion 
given by Eq. (1) and (2). In Section III, we reformulate the R-matrix theory 
based on the Hilbert-Schmidt theorem to obtain the R-function in terms of 
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a series which is uniformly and absolutely convergent. In Section IV, we 
discuss methods of improving the convergence of the R-function series when 
it is truncated for practical applications. A summary and conclusions are 
given in Section V. 

II. R-function 

In order to describe the formal procedure employed in the R-matrix the- 
ory, we consider the simplest case of potential scattering for spinless particles 
with only the elastic scattering channel being open. 

The radial wave function ue{r) in the interior region < r < a satisfies 
the Schrodinger equation 



d^ue ^ 2ij,V{r) ^ i{i + 1) 



+ (^T2^ + -^^)Mr) = k'uiir), (3) 



where /i is the reduced mass, V{r) is the interaction potential in the internal 
region < r < a, and k"^ = 2n,E/fi^. 

In the conventional theory [3, 18, 37], Ui{r) in the internal region (0 < 
r < a) is expanded in terms of complete set of states «^(r) given within the 
region < r < a. These states are the solutions of the equation 

' + (-^^ + = A;M(r), (4) 



satisfying the R-matrix boundary conditions 

uiiO) = 0, 



a I dul^ 

' )r=a 



(5) 



li^(a) dr 
and the orthonormality conditions 

r ui{r)ui,{r)dr^6,y. (6) 

Jo 

In the region < r < a, U£{r) may be expanded in terms of the eigen- 
functions Ux{r). 

oo 

Mr) = E ^Wxir), (0 < r < a), (7) 

A=l 
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where 

drue{r)ux{r). (8) 

As we show below, cither (i) the expansion (7) does not converge uni- 
formly, or (ii) term-by-term differentiation is not admissible [3, 18], or both 
(i) and (ii) may be applicable. Prom Green's theorem [3, 18] and the bound- 
ary conditions (5), we find 

f 1 uUa) , due „ , 

Substitution of Eq. (7) into Eq. (6) gives 

«,(r)[a^ - Bu,];l^ = - E (10) 



If we now define 



and assume that 



ui{r)ui{a) _ - {ui{a)r .... 

A=l '^A A=l ''^A 

we find that R^^^ relates the amplitude of U£ to its derivative at the boundary 
by the relation 

R^^^ ^ ue{a)[a^ - Bue];l^. (13) 

Once i?^^) has been calculated, the K-matrix and cross-section can be easily 
determined. 



From the fact that 



and 



A=l '^^A ^ 



nr^ „A _/ ue{r)Yi{r'), r < r', . 



4 



with 

where y(,{r) is the irregular sohition of Eq. (3) with boundary conditions 
= 0, and u'f^Uf^ — y[uf = 1, it can be seen that Eqs. (7) and (10) can 
be obtained from the spectral decomposition, Eq. (14a), of Green's function 
Ge{r,r'). This bilinear series, Eq. (14a), converges in L2(0 < r < a). 

We note that the completeness of the states does not guarantee vahdity 
of Eq. (12). To demonstrate this statement, let us first consider a special 
case of boundary conditions S = oo, or 

uiiO) ^ uiia) ^ 0. (16) 

In this case, we can find from Eq. (9) that 

J _ u,{a) dui 

- ^23^ ^^^'■=«' ^^^^ 

and substitution of Eq. (17) into Eq. (7) gives 

Mr) _ ^ ujir) , du{ 

If one tries to obtain this value of hm M^) ^ from the right side of Eq. 

ue{a) 

(18) taking the limit term by term, one obtains a null result, because of Eq. 
(16). 

In the case of the boundary conditions (5) we can obtain from Eq. (10), 
that 

where 

V{r) = [«-J^ - Bue{r)]/[a^ - Bu,]r=a. (20) 

Once again wc obtain a null result for ri{a) = 1 by differentiating term by 
term and taking the limit term by term of the sum in Eq. (19) and using Eq. 
(5). The explanation for these paradoxes is that either (i) the expansion (7) 
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or (ii) its derivative series, obtained by differentiating tlie individual terms 
of the expansion (7), is not uniformly convergent in the neighborhood of the 
surface. Or they may be due to both (i) and (ii). This difficulty associated 
with the expansion (7) has been known for many years [3,18]. 

III. The Hilbert-Schmidt theorem formulation of the R-matrix the- 
ory. 

Let us rewrite Eq. (3) in an integral form 



Ui{r) = Mr) + {k' - / Ke{r, r')ue{r')dr' , (21) 

Jo 

with 

KArr')--l ^'^'^^'^''^^ " ^ (22) 
\Ye{r)X,{r'), r'<r, ^^^^ 

where Xg{r) and Yg{r) are regular and irregular solutions, respectively, of the 
following equation 

and satisfy the following conditions 

X,(0) = 0, 



Ye{a) dr 

and 



a dYe _ 



(24) 



^x.(r)-^r.(0^-i. 

dr dr 

B is the same as one the introduced in Eq. (5). is an energy independent 
constant satisfying a condition 

«V^A (A = 1,2, ....), (25) 

and 4>i{r) is related to X^lr) by 

Mr)^aXe{r), (26a) 
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where a is an energy dependent constant given by 

a = [(a^ - Bu,{r)]r=a/[a^^^ - i?X,(r)].=„. (266) 

The integral equation (21), which is not the Lippmann-Schwinger type equa- 
tion, was first introduced in [9] for the — 0, B — cx) case. Eq. (21) has a 
unique solution, since 



pa pa 

/ / Kl{r, r')dr dr' <oo, 
Jo Jo 



i.e. Ki{r, r') is completely continuous and self-adjoint kernel [39]. Let jx{X — 
1,2,...) be eigenvalues of the Hermitian continuous kernel K£{r, r') 

<(0=7a r K,{r,r')ui{r')dr', (27a) 

Jo 

with 

lx = kl- k\ (276) 

As it is well known, the eigenvalues 7a are real, and the functions Ui{r) and 
M^(r) are continuous. Due to the Hilbert-Schmidt theorem [39], the following 
expansion 

/ Ury)u,{r')dr' ^ ^ cWx{r) (28) 

•'° A=l 

converges uniformly and absolutely over < r < a, and, if k"^ ^ k^, the 
unique solution Ui{r) of the integral equation (21) appears in the following 
form of a series which is uniformly and absolutely convergent over <r < a 
(by Schmidt's formula): 

Mf.^-R,/ 1-1 - ^^(^) k^-n^ ^ ui{r)ui{a) 

ue[r)[a^^ ""^^^^ " [a^M _ 5x,(r)] « kikl-^^m-k^Y 

(29) 

If we now define 

R^^\k') = R^^\^') + ^-^^ t (30) 



7 



where 

R^^\k')= , (31) 

we find that R^^\k'^) relates the amplitudes ui to its derivative on the bound- 
ary by the relation (13). Because the series (29) converges uniformly and 
absolutely, the following equation 

is vahd [40], and hence the expansion (29) is free of difficulties encountered in 
the expansions given by Eqs. (7) and (10). Series (30) can be also obtained 
from the dispersion formula (11) by separating the energy independent term 
R^^\k^) in the R-function, Eq. (11), with satisfying condition (25). Our 
derivation has shown that the dispersion expansion (30) converges absolutely, 
and exhibits the general energy-dependence of the R-function. The expansion 
given by Eq. (29) is a main result of this paper. The proof of the absolute 
convergence of the series (30) in case = was given by M. Schiffer and V. 
Bargmann. Their proof is reproduced in [41]. 

IV. Improving the Convergence 

In general, the R-function has an infinite number of pole terms. According 
to Courant's minimax considerations, if V{r) is bounded, no kj^ differs from 

the corresponding value of k^, {k^^^^, for noninteracting case {V{r) = 0) by 
more than the bound [3]. Consequently, the general term of the series (11) 

for fixed /c^ behaves as 1/A^ since {kf^^ oc A^, while the general term of the 
series (30) behaves as 1/A^ as A — > oo. For the case of the dispersion formula 
(11), truncation of the R-function by a finite number (N) of terms gives 

While for the case of the dispersion formula (30), we have 

<'(*^) = e' + i|;gi^, (34) 

where 

i;r=fl(«'M-i|:|f^. (35) 
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The general method of improving the convergence is to separate and sum 
the slowly converging parts of the series [42] . It is obvious that there are many 
possibilities to obtain a rapid convergence. For example, the expansion (30) 
can be represented in the form 

[Kq — K ) 

(36) 

where k^ is an energy independent constant [k^ ^ n'^jk^ ^k\^\ = 1,2,...). 
Expansion (36) converges much faster than (30) (general term behaves as 
1/A^), and truncation of Eq. (36) by a finite number of terms gives 

ieV-(^^)^sr+;^^«r'+is^. (37) 



where 



and 



(38) 



However, for the case of Eq. (37), we have introduced an additional param- 
eter Ri^\ and we do not expect a weak dependence of Eq. (37) on this 
parameter r[^\ 

To obtain a faster convergence, we introduce a trial potential V{r) and 
remove the corresponding R-function R^^\k'^) obtained with V{r): 

R^^\k-) = R^^\n-) + ^ t ^""llf^ (40) 

« t^i {kl-K^){kl-k^) 
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For the case of the dispersion formula (11), this method has been used in 
many papers [11, 15, 18]. For the case of Eq. (30), we have 

^ e-K' g , {ui{a)r {ui{a)f ^ 



(41) 

It can be shown (see Appendix) that the general term of Eq. (41) behaves 
as 1/A^ for any bound V{r), and hence we expect a weak V{r) dependence 
for the following approximation 

^iji"' + i £ ^1^. (42) 

A=l 



where 



g,{e) = R^-\e) - - - E (43) 

^ A=l '^X~ ^ 



^r=i?(-)(.^)-l E (44) 

^ A=l '''A 

and ui{r) are solutions of Eq. (4) with the trial potential V{r). Note that 
the case of V^(r) = was considered in [9, 43]. For practical calculations 
of the R-matrix with the approximation (42) for an incident nucleon, it is 
possible to use a simple squared-well potential 

V{r)^-Vo9{a-r), (45) 

where Vq = . The wave number Kq is independent of the mass number 
A and is approximately the same for all nuclei {Kq ^ Ifm) [44]. 

V. Summary and Conclusions 

Using the Hilbert-Schmidt theorem and the integral equation, Eq. (21), 
wc have reformulated the R-function theory in terms of the expansion given 
by Eq. (29) which is uniformly and absolutely convergent for all values of 
< r < a. This expansion, Eq. (29), can be differentiated term by term 
in the neighborhood of the surface. Our reformulation solves the existing 
formal problem of how to formulate the R-matrix theory without the use 
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of expansions which are not uniformly convergent. A possible method for 
improving the convergence of the R-function series is given when the series 
is truncated for practical applications. 
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Appendix 



In this Appendix, we show that the general term of Eq (41) behaves as 
1/A^. We rewrite Eq. (4) in a form of the Volterra integral equation 



u\[r = 0n^[r 



) + r K(r,r')u{(r')dr', (A.l) 
Jo 

K{r, r') = [x^oA(r)/oA(0 - ct>Ur)xUr')\^i.r% (A.2) 



where 
with 



(f)^^ in Eq. (A.l) is a regular solution of the Schrodinger equation for 
noninteracting case 



satisfying the R-matrix boundary conditions (5) and the orthonormality con- 
ditions (6). XoA ™ Eq. (A. 2) is the irregular solution of Eq. (A. 3) 



XoA(r) = 0'oA(r) [[<PU^)]-'dx, 

Jo 



and k1 is defined from the condition 



r Xox{rMr)ui{r)dr = 0. {AA) 
Jo 



For any bound and continuous V{r), K{r, r') is also continuous and bound, 
and hence the Neumann series (iteration series) 

CXD 

p=0 

converges uniformly and absolutely over < r < a [39], where is a 
product of the operators K and the function K{r, r') is the kernel of the 
linear integral operator K. 
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^From Eq. (A. 5), we can obtain 

lim «i(a)=0^oA(«) + O(^), (A6) 

A— *oo A 

and hence the general term of Eq. (41) behaves as 1/A^ for any bound and 
continuous V{r) and V{r). 

We note an important fact that wc do not need "smallness" of V(r) for the 
convergence of the Neumann series, Eq. (A. 5), in contrast to the conventional 
perturbation expansion. 
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